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We present a guiding principle for designing fermionic Hamiltonians and quantum Monte Carlo 
(QMC) methods that are free from the infamous sign problem by exploiting the Lie groups and Lie 
algebras that appear naturally in the Monte Carlo weight of fermionic QMC simulations. Specifi¬ 
cally, rigorous mathematical constraints on the determinants involving matrices that lie in the split 
orthogonal group provide a guideline for sign-free simulations of fermionic models on bipartite lat¬ 
tices. This guiding principle not only unifies the recent solutions of the sign problem based on the 
continuous-time quantum Monte Carlo methods and the Majorana representation, but also suggests 
new efficient algorithms to simulate physical systems that were previously prohibitive because of 
the sign problem. 
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One of the biggest challenges to classical simulation of 
quantum systems is the infamous fermion sign problem 
of quantum Monte Carlo (QMC) simulations. It appears 
when the weights of configurations in a QMC simulation 
may become negative and therefore cannot be directly 
interpreted as probabilities [1], In the presence of a sign 
problem, the simulation effort typically grows exponen¬ 
tially with system size and inverse temperature. 

While the sign problem is nondeterministic polynomial 
(NP) hard [2], implying that there is little hope of finding 
a generic solution, this does not exclude ad hoc solutions 
to the sign problem for specific models. For example, one 
can sometimes exploit symmetries to design appropriate 
sign-problem-free QMC algorithms for a restricted class 
of models [3]. However, it is unclear how broad these 
classes are and it is in general hard to foresee whether a 
given physical model would have a sign problem in any 
QMC simulations. The situation is not dissimilar to the 
study of many intriguing problems in the NP complexity 
class, where a seemingly infeasible problem might turn 
out to have a polynomial-time solution surprisingly [4], 

A fruitful approach in pursuing such specific solutions 
is to design Hamiltonians that capture the right low 
energy physics and allow sign-problem-free QMC sim¬ 
ulations at the same time, called “designer” Hamiltoni¬ 
ans [5]. This naturally calls for design principles. For 
bosonic and quantum spin systems a valuable guiding 
principle is the Marshall sign rule [6, 7] which ensures 
nonnegative weight for all configurations. The design of 
the sign-problem-free fermionic Hamiltonians is harder. 
The method of choice for fermionic QMC simulations 
are the determinantal QMC approaches, including tra¬ 
ditional discrete-time [8] and new continuous-time ap¬ 
proaches [9-13]. Both approaches map the original inter¬ 
acting system to free fermions with an imaginary-time 
dependent Hamiltonian. The partition function is then 
written as a weighted sum of matrix determinants after 


tracing out the fermions [8, 9, 12]: 

Z = J2fc det [/ + Te~ H drffc(r) 


(1) 


where fc is a c-number and Hq(t ) is an imaginary-time 
dependent single-particle Hamiltonian matrix (whose 
matrix elements denote hopping amplitudes and onsite 
energies on a lattice), both depending on the Monte Carlo 
configuration C. T denotes the time-ordering and / is the 
identity matrix. The appearance of the matrix determi¬ 
nant complicates the analysis of the sign problem because 
it is often not straightforward to see the sign of the Monte 
Carlo weight of a given configuration [14, 15], and the 
sign of the determinant is related [16] to the Aharonov- 
Anandan phase [17] of the imaginary-time evolution. The 
situation is further complicated by the fact that even for 
a given physical model the choice of the effective Hamil¬ 
tonian He is not unique (it depends on details of the 
QMC algorithm such as whether and how to perform an 
auxiliary field decomposition) and the specific choice may 
affect the appearance of the sign problem [14, 18, 19]. 

One successful guiding principle for fermionic simula¬ 
tions that has been discovered in the context of nuclear 
physics [20, 21], lattice QCD [22] and condensed matter 
physics [23] relies on the time-reversal-symmetry (TRS) 
of the effective Hamiltonian He . TRS ensures a nonneg¬ 
ative matrix determinant in Eq. (1) because the eigen¬ 
values of the matrix necessarily appear in Kramers pairs. 
A typical example of this kind is the attractive Hubbard 
model at balanced filling of two spin species, where after 
decomposition of the interaction term the Monte Carlo 
weight even factorizes into product of two identical ma¬ 
trix determinants. Additional conditions such as half Fil¬ 
ing and bipartiteness of the lattice lead to a solution of 
the sign problem for the repulsive Hubbard model. See 
Refs. [23, 24] for a thorough discussion and Refs. [25-27] 
for several recent applications of the TRS principle. 

Unfortunately, besides the quite intuitive TRS princi- 
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pie [21-23], a broad criterion for the sign of the matrix 
determinant is still lacking. Recent progresses on solving 
the sign problem in a class of fermionic models using the 
continuous-time quantum Monte Carlo approach [28] and 
the Majornara representation [29] provide hints about 
such a guiding principle. For example, one could search 
for real-antisymmetric matrices with nonnegative deter¬ 
minant [28, 30], or try to split the fermionic operator 
into Majorana fermions for a potential cancellation of the 
sign [29]. However, compared to the TRS principle [21- 
23] , both approaches are still not enlightening enough to 
serve as a guiding principle. Moreover, because of the dif¬ 
ferent appearances of the two solutions [28, 29], it is un¬ 
clear what is the connection between them and whether 
there is a deeper underlying reason for such solutions. 

In this Letter, we present a guiding principle that not 
only unifies the two recent solutions to the sign prob¬ 
lem [28, 29], but also suggests a general strategy that 
enables us to discover solutions to the sign problem for a 
broader class of fermionic models. The guiding principle 
exploits the symmetry of the effective Hamiltonian He 
and consequently the Lie group structure of the evolu¬ 
tion matrix Te~ h drH c j n particular, the split orthog¬ 
onal group 0(n, n ) is formed by all 2 n x 2n real matrices 
that preserve the metric r] = diag(l,..., 1 , — 1 ,..., — 1 ) 

n n 

M T r)M = T]. (2) 

Similar to the Lorentz group 0(3,1), a familiar example 
in relativistic physics, the 0(n, n ) group contains four 
components. More explicitly, writing the matrix M in 

the form M = ( ^ 12 ) with n x n sub-blocks, one 

\ M 21 M 22 J 

has |det(Mn)| > 1 and |det(M 22 )| > 1 [31]. The four 
components of 0(n, n ) can be classified by the signs of 
det(Mn) and det(M 22 ), denoted as 0 ±:t (n,n). Different 
components can only be connected by improper rotations 
that change the sign of the determinant of the sub-block 
Mu or M 22 • Only the 0 ++ (n, n) component forms a 
subgroup because it contains the identity element. 

Theorem If M belongs to the split orthogonal group 
0(n,n), then the following statements hold [32-34] 


det (/ + M) 


> 0 , 

if M e 0 ++ (n,n), 

(3a) 

< 0 , 

if M € O (■ n,n), 

(3b) 

= 0 , 

otherwise. 

(3c) 


This rather strong statement about the definite sign of 
the matrix determinant, no matter positive or negative, is 
invaluable for the determinantal QMC simulations. Fur¬ 
thermore, we have the following 

Corollary Given an arbitrary number of real matrices 


Ai that satisfy r/Air) = —Aj, we have [32, 33] 



( 4 ) 


The proof follows immediately by noticing that Ai lies 
in the Lie algebra of the group 0{n,n ) [35]. Each fac¬ 
tor of the matrix product e Ai is an exponential from 
the Lie algebra to the Lie group 0(n, n ) in the identity 
component, thus Eq. (4) is a consequence of Eq. (3a). 
Note that the form of the matrix determinant of Eq. (4) 
resembles the weight that appears in the determinantal 
QMC calculations Eq. ( 1 ) [ 8 , 9, 12], 

Before moving on, we comment on the general rele¬ 
vance of Eq. (3) and Eq. (4) to physical problems. On a 
bipartite lattice, the parities of the sublattices naturally 
provide the metric 77 appearing in Eq. (2). To further 
reveal its physical meaning, we write an element in the 


Lie algebra Ai = 


explicitly. In the special 


Ci Bi 

Bf A 

case ol Ci = Di = 0, Ai can be recognized as a bi¬ 
partite single-particle Hamiltonian and the condition on 
A, has appeared in Eq. (4) of Ref. [28]. The Corollary 
(4) states that the partition function of such bipartite 
imaginary-time dependent noninteracting system is non¬ 
negative [36]. Moreover, in general the matrix Ai does 
not need to be symmetric. The condition on Ai only re¬ 
quires Cf = —Ci and Df = — Di , thus provides more 
flexibilities in designing the QMC approaches. 

To see how the above rigorous mathematical state¬ 
ments apply to determinantal QMC simulations of phys¬ 
ical systems, we consider first the spinless t — V model 
on a bipartite lattice, 


H — 'y ' c\Kijdj 


(hi) L 


V 


hi + hj 




-r 


■ ( 5 ) 


Here c\ and Ci are fermion creation and annihilation op¬ 
erators and hi = ejei is the occupation number operator 
on site i. There are 2 n lattice sites which split into two 
sublattices A and B. In accordance with the metric rj , we 
sort the sites by placing all sites in A before those in B. 
The bipartite hopping matrix K has zeros on the diagonal 
and is real-symmetric, therefore fulfills the requirement 
r/Kr] = — K t of the Corollary [28]. The second term 
is a repulsive interaction between nearest-neighbors (i,j) 
(belonging to different sublattices) and we introduced a 
constant shift T, which will play a crucial role in later 
discussions. 

We employ the continuous-time quantum Monte Carlo 
(CT-QMC) framework [9-12] in the following analysis. 
This approach is free from time discretization errors, 
as efficient [12] and more flexible and powerful [13, 37] 
than the discrete-time counterpart [ 8 ]. Furthermore the 
discrete-time algorithms can be derived as a restricted 
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version of CT-QMC on an equidistant grid of imaginary- and perform an expansion in the interaction term [9] 
times [38] and our results apply to them as well. We 
rewrite the Hamiltonian Eq. (5) as H = H 0 + Yl(ij) ®ij 

I 


OO 

(^")=E £ ■ 

■■ E J dn.. 

f0 

■ / dr k Tr 

e (/? fk)Ho (~Vikj k ) •• 

■(— Kh) e TlH ° 

(6) 

fc=o (uji) 

(ikJk) ° 

J Tk- 1 





At this point there are multiple ways to proceed, which 
result in distinct CT-QMC algorithms, differing in ef¬ 
ficiency and the use of auxiliary fields (see [13] for an 
overview). In particular, we choose the following aux¬ 
iliary field decomposition of the interaction term to re¬ 
veal the connections of various solutions to the sign prob¬ 
lem [28, 29], 


~Vij = 77 E exp 

<T=± 


a\ (c\dj + 


(7) 


where A = acosh (l + jp) is a real number for repulsive 
interaction V and positive shift T. The decomposition 
Eq. (7) is valid because the operator <5 = (c\cj + cjc,) 
satisfies o 2 = o 4 = n* +nj — 2fiihj when i ^ j. Compared 

I 


OO 


Z = E 

/c—0 



£ £-£ jfV- 

1=± o-fe=± 


where the matrix (A fj)i m = aX(8u5 m j + SijS m i) accord¬ 
ing to the exponential factor of Eq. (7). Equation (8) 
is in the general form of Eq. (1) and the matrix deter¬ 
minant has the form of Eq. (4). The interaction vertex 

r , ... ( cosh A usinhA 

e ^ performs a hyperbolic rotation I ypjj y cosh A 

in the relevant 2x2 block involving the sites i,j. Im¬ 
portantly, both the original hopping matrix K and the 
auxiliary Hamiltonian matrix AA satisfy the condition of 
the Corollary (4). The weight (8) is therefore nonneg¬ 
ative and there is no sign problem. The Monte Carlo 
method can be used to sample the summations over the 
interaction bonds and the auxiliary fields as well as the 
integrations over the imaginary times on equal footing, 
see Refs. [12, 13] for details about efficient Monte Carlo 
simulation of Eq. (8). 


Using the auxiliary field to decouple the interaction 
vertex is not the only way to formulate a sign-problem- 
free QMC approach for the model (5). The Theo¬ 
rem (3b,3c) apply to other components of the 0(n, n) 
group and connect the above solution to the solu¬ 
tions based on the continuous-time interaction expansion 


to the conventional decompositions routinely employed in 
the determinantal QMC simulations [9, 39], the auxiliary 
field in Eq. (7) couples to fermion hoppings instead of 
the density operators [40, 41] . This is one of the key 
ingredients to avoiding the sign problem. In retrospect, 
this choice can be motivated by the Corollary (4). 

Plugging Eq. (7) into Eq. (6), the square bracket be¬ 
comes a product of exponentials of fermion bilinear oper¬ 
ators. The trace therefore acquires an appealing physical 
meaning: it is the partition function of an imaginary-time 
dependent noninteracting system, which evolves alterna¬ 
tively under the free part of the original Hamiltonian Hq 
and hopping with an amplitude aX between the sites z, j 
that belong to different sublattices. Tracing out these 
free fermions, one obtains 


/A 

/ dr k det 

I + e -^~ Tk)K e A S'* . 

..e A Sie- riK 

' Tk- 1 




( 8 ) 


method (CT-INT) [28, 30] and the related but more ef¬ 
ficient LCT-INT method [12, 13]. These methods corre¬ 
spond to special choices of the shift T = —V/A [42] which 
results in a purely imaginary coupling strength A = zzr 
in Eq. (7). The vertex matrix e A o thus has the form 

^ in the relevant 2x2 block independent of 

the auxiliary field, which is equivalent to rewriting the 
interaction term vij = h e »(ni+"j) i n the LCT-INT ap¬ 
proach [12, 13]. The vertex matrix maps the evolution 
matrix in Eq. (8) back and forth between the 0 ++ (n,n) 
and the O (n,n) components, because the sites i,j be¬ 
long to different sublattices and the vertex matrix flips 
the signs of both det(Mn) and det(M 22 )- The matrix 
determinant in Eq. (8) is thus non-positive for an odd 
number of vertices according to Eq. (3b). However, the 
negative value T = — V/A cancels this sign due to a pref¬ 
actor (^) fc in the weight. Hence Theorem (3) ensures the 
absence of a sign problem [43] in the auxiliary-field-free 
(L)CT-INT simulations [13, 28, 30]. 

To take full advantage of the Corollary (4), one 
can further consider long-range interactions in the 
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model Eq. (5), e.g., attractive interaction between sites 
belonging to the same sublattice [28, 29]. We decouple 
the interactions as 


~Vij = 77 exp aX ~ £ l £l ) 


cr—=b 


(9) 


The coupling strength A = acos (l + is real for at¬ 
tractive interactions and any positive shift T > |U|/4. 
The effective single-particle Hamiltonian in the exponen¬ 
tial of Eq. (9) is antisymmetric and connects sites in 
the same sublattice, thus fulfills the requirement of the 
Corollary (4). There is no sign problem either [44]. 

Alternatively, in Eq. (7) and Eq. (9) one can split 
a fermion into two Majorana operators [29] and iden¬ 
tify two complex-conjugate factors in the Monte Carlo 
weight [45]. It is however clear that the unconventional 
decoupling in the hopping channels in Eqs. (7,9) to re¬ 
spect the Corollary is the underlying reason of a nonneg¬ 
ative matrix determinant. In light of Eq. (4), rewriting 
the fermions using Majorana operators is unnecessary in 
the Monte Carlo simulations. Nevertheless, the Majo¬ 
rana representation [29] is an ingenious way to prove the 
Corollary. 

We have shown that Theorem (3) unifies the recent so¬ 
lutions of the sign problem [28, 29] as different choices of 
the constant shift T. The Corollary (4) is particularly in¬ 
structive as it suggests that one just needs to decompose 
the original interacting model into free effective Hamil¬ 
tonians that satisfy the condition of Eq. (4) in order to 
avoid the sign problem. The mechanism of solving the 
sign problem using Eq. (3) and Eq. (4) goes beyond the 
previous understandings based on the TRS principle [21- 
23]. This can be easily seen from the fact that the real 
eigenvalues of the matrix I + M are not necessarily dou¬ 
bly degenerate as required by the Kramers theorem [46] . 

As a further application [32] we consider the following 
two-flavor Hubbard model on a bipartite lattice 

H = J2 Yu^ K r> £ 3« + Hw 

a={T,U i,j i 

Vi = U (nifhii - _ r, (10) 

where the real hopping matrix K a connects the same fla¬ 
vor a on different sublattices. The model Eq. (10) covers 
a variety of interesting physical systems that were pre¬ 
viously inaccessible for determinantal QMC simulations. 
For example, the choice = rK "1 with a ratio 0 < r < 1 

realizes the asymmetric Hubbard model, which was im¬ 
plemented recently in a one-dimensional optical lattice 
with a tunable ratio r [47], On the other hand, one can 
also choose to have spatially anisotropic hopping ampli¬ 
tudes for each flavor, therefore to realize Hubbard models 
with mismatched Fermi surfaces [48]. 

All these cases break the SU(2) spin symmetry as well 
as the time-reversal symmetry, therefore are not guar¬ 


anteed to be sign-problem-free according to the conven¬ 
tional TRS principle [21-23]. However, one can now 
solve the sign problem using the insights provided by 
the Corollary (4). We first consider the U > 0 case for 
simplicity. Enlightened by the new understanding, we de¬ 
couple the interaction term Eq. (10) similarly to Eq. (7) 
and obtain an auxiliary field coupled to the local spin-flip 
(cj^cq, + c^Ci-)-), which connects different flavors on the 
same site. Thus for the ordering of the spin-orbital (A t, 
B B A \), it is easy to see the effective Hamilto¬ 
nians are bipartite and symmetric, therefore satisfy the 
condition of the Corollary. This shows that an auxiliary 
field coupled to the x-component of the spin operator is 
sign-problem free for the model (10) [49]. The attractive 
case can be studied without a sign problem by perform¬ 
ing a particle-hole transformation to the model. Alter¬ 
natively, one can perforin the decomposition according 
to Eq. (9) for attractive interactions, thus have a sign- 
problem-free simulation with the auxiliary field coupled 
to the y-component of the spin operator. Moreover, there 
is no sign problem even when we explicitly add spin-flip 
terms in the Hamiltonian as long as the hopping matrix 
satisfies the condition of Eq. (4). This covers a large 
class of compass Hubbard models [50] which are relevant 
to multiorbital and ultracold atoms systems [51-53]. 

Using the special choice of F = —C//4, the above so¬ 
lution reduces to the (L)CT-INT formulation and the 
determinant of the two flavors factorizes into two parts 
in the absence of the single-particle spin-flip terms. Even 
though the two determinants are not necessarily equal 
due to the broken TRS, Theorem (3) ensures that they 
have the same sign because the evolution matrix of the 
two flavors lie in the same component of 0(n, n). In con¬ 
trast to the case of spinless fermions, the vertex matrix of 
Vi = lL e ™(n i1; +n il ) can bring the evolution matrix into all 
four components of the 0(n, n) group since each vertex 
matrix changes the sign of either det(Mn) or det(M 22 ) 
of both flavors. The Monte Carlo weights of odd ex¬ 
pansion orders vanish because of Eq. (3c). Although the 
matrix size in the LCT-INT simulation is only half of the 
previously discussed auxiliary field approach, the use of 
two-vertices insertion/removal updates [10] in the Monte 
Carlo simulation leads to more complicated updates and 
measurement procedures [54, 55] . The auxiliary field ap¬ 
proach may thus be advantageous. 

These solutions to the sign problem can also be applied 
to projector QMC methods [13, 56, 57] which sample 
the ground state wave-function overlap {'t , T\e~ eH I^t) 
instead of the partition function. One can choose the 
trial-wave function |4 /t) as the ground state of a single¬ 
particle trial Hamiltonian that fulfills the condition of 
Corollary (4) to avoid the sign problem. 

All the sign-problem-free models solved by Eq. (3) and 
Eq. (4) in this Letter are at half-filling on bipartite lat¬ 
tices with particle number conservation [58]. It will be 
interesting to see whether one can even go beyond this 
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constraint. Conversely, we emphasize that the require¬ 
ments of Eq. (3) and Eq. (4) are by no means the neces¬ 
sary conditions for a sign-problem-free QMC simulation. 
There should be more “de-sign” principles of this kind 
for fermionic Hamiltonians and quantum Monte Carlo 
methods. Our work suggests it is fruitful to exploit the 
inherent Lie group and Lie algebra structure in the Monte 
Carlo weight to search for such “de-sign” principles. In¬ 
cidentally, both the split orthogonal group and the TRS 
“de-sign” principle seem to be related to the ten-fold way 
classification of random matrices [59] . It would be inter¬ 
esting to generalize them to other symmetry classes [60- 
62] and draw connections to recent topological classifica¬ 
tion of gapped free-fermion systems [63-65] . 

Furthermore, findings reported in this paper apply as 
well to fermions coupled to quantum spins or Z 2 gauge 
fields. The Theorem (3) ensures a matrix determinant 
with a definite sign after integrating out fermions as long 
as the split orthogonal group structure is respected. This 
allows to design new sign-free models relevant to lattice 
gauge theories [32]. 
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Proof of Eq. (3): A showcase of the open science 
movement 

The Theorem Eq. (3) was proved by Gergely Harcos 
and Terence Tao, which was triggered by Lei Wang’s con¬ 
jecture posted on MathOverflow. Below we review the 
history of the proof. 

On May 1st 2015, Wang posted a weaker version of 
the Corollary Eq. (4) on MathOverflow formulated as 
the following 

Conjecture Given an arbitrary number of real matrices 
A i = ( _B t "o* ) ’ det + n, e Ai ) holds. 

Tao soon pointed out the matrix product JX, c lies 
in the split orthogonal group. Christian Rending and 
Will Sawin confirmed this observation and clarified the 
Lie algebra of the 0(n,n ) group. A “counterexample” 
to the conjecture (a matrix M £ 0(2,2) which has 
det(/ + M) < 0) was posted but realized later to lie in 
the O (2, 2) component. These discussions helped to 
clarify the topological structure of the split orthogonal 
group. Harcos then formulated the question in the form 
of Eq. (3a) and outlined a strategy of the proof. Tao 
finalized the proof and summarized a streamlined ver¬ 
sion in his blog (Proposition 3 in [34]). Wang suggested 
Eq. (3b) and Eq. (3c) as further extensions in the com¬ 
ment section of Tao’s blog and got affirmative answers 
from Tao and Harcos. There are other people who con¬ 
tributed to the discussion and confirmation of the proof, 
some are anonymous users of MathOverflow. 

Three-flavor Hubbard model 

Consider the following three-flavor Hubbard model [28] 

11 = E E 4^-a +J2J2 * 

a={l,2,3} i,j i a<0 

~a/3 _ ,r (~ « hi a + flip \ 

I j — U a/3 l riiaTlip ^ 1 fa/3* (ff) 


Again, the hopping amplitude K a is real and bipartite, 
but it does not need to be the same for each flavor. With 
attractive onsite interactions U a p < 0 a decomposition of 
the interaction term similar to Eq. (9) gives an auxiliary 
field that couples to (c\ a Cip — c\pCi a ). For the ordering of 
the spin-orbital (Al, *42, A3; HI, £>2, £>3) the exponential 
factor forms a real antisymmetric matrix that fulfills the 
condition of the Corollary. There is no sign problem in 
the simulation. Alternatively, similar to the discussions 
in the main texts, the determinant factorizes for the three 
flavors with the choice of T a p = — U a p/A. Their prod¬ 
uct is nonnegative and the formalism reduces the sign- 
problem free (L)CT-INT simulation [28]. 


Z 2 Lattice gauge theory 

Consider Z 2 lattice gauge theory [66] coupled to N- 
flavors of fermions 


( i,j) plaquettes 

N 

( 12 ) 

(i,j) a=l 


where a x and a z are Pauli matrices defined on links of a 
d-dimensional bipartite lattice. The fermions are at half¬ 
filling and live on lattice sites. The first line of Eq. (12) 
describes the Z 2 lattice gauge theory in d+1 dimensions. 
The second line couples the Z 2 gauge fields to fermion 
hoppings across the links. The model is sign problem 
free because after integrating out the fermions the ma¬ 
trix determinants satisfy the Corollary (4). Notice that 
the Monte Carlo weight is nonnegative for any flavor N 
because one does not rely on cancellation of signs between 
different flavors. 



